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摘要 :我 国 西南 喀斯特 地 区 石 漠 化 面积 已 实现 净 减 少 ,植被 状况 具有 明显 改善 。 为 轩 更 清晰 的 了 解 该 区 植被 变化 情况 及 其 影响 
因素 的 区 域 差 异 ,采用 长 时 间 序 列 遥 感 数据 ,综合 运用 空间 自 相 关 分 析 . 主 成 戎 分 析 (BECA ) 和 地 理 加 权 回 归 ( GWR ) 等 研究 方 
法 ,分 析 生 态 工程 实施 以 来 汗 桂 黔 喀斯特 植被 变化 及 其 主要 影响 因素 的 空间 非 平 稳 性 。 结 果 表 明 : 与 1982 一 2000 年 相 比 ， 
2001 一 2011 年 生长 季 归 一 化 植被 指数 (CSN ) 在 整个 研究 区 域 都 有 增加 且 具 有 显著 的 空间 集聚 性 ( 到 为 0.90) ,但 增加 程度 在 空 
间 上 具有 差异 性 ( 变异 系数 为 43%) ;影响 汗 桂 黔 植被 变化 的 主要 因 奈 甸 括 氧 候 因子 .土壤 质地 、 人 类 活动 ,水 分 有 效 性 .土壤 养 
分 和 社会 经 济 条 件 , 且 对 植被 变化 的 影响 程度 随地 理 位 置 的 变化 而 变化 。 不 同 工 程 地 貌 类 型 区 内 ,影响 植被 变化 的 主导 因素 不 
同 , 且 存 在 显著 的 空间 差异 性 , 需 综合 考虑 植被 变化 主导 因素 的 区 域 差异 来 调整 或 改进 后 续 生 态 工程 措施 。 

关键 词 : 汗 桂 黑 ;植被 变化 ;地 理 加 权 回 归 ; 空 间 非 平稳 性 ; 主导 因素 
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Abstract: With the implementation of ecological restoration, the net area of rocky desertification has been reduced and 
vegetation coverage significantly improved in karst areas, southwest China. It is crucial to explore the trends in vegetation 
change’ and distinguish the driving forces under different karst terrain conditions. Based on the long-term time-series of 
GIMMS NDVI 3g images, the present study aimed to identify the vegetation changes at different spatial and temporal scales 


and non-stationary relationships between vegetation changes and driving forces in karst regions of Yunnan, Guangxi, and 
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Guizhou Provinces with a combination of spatial autocorrelation, principal component analysis (PCA), and geographically 
weighted regression. The results showed that from 2001 to 2011 the growing season normalized difference vegetation index 
(GSN) significantly increased compared to that from 1982 to 2000 in the three regions. The GSN appeared spatially 
assembled (the Moran's index J, was 0.90), while the extent of increase was significantly diverse for different areas (the 
variable coefficient was 43%). The PCA indicated that the main driving forces of vegetation change were climatic factors, 
soi] texture, human activities, water availability, soil nutrients, and the level of social and economic conditions. The extent 
of impacts of the driving forces on vegetation change varied with different geographic positions. In addition, the mainedriving 
forces were spatially different under different karst landforms. Our study indicated that the differences of vegetation drivihg 
forces should be taken into consideration during future implementation of ecological restoration under different karst 


landforms. 
Key Words: karst region; vegetation dynamics; GWR model; spatial non-stationary; dominant factors 


为 遏制 石 漠 化 土地 扩展 、 加 快 石 漠 化 土地 治理 进程 ,2008 年 ,国务 院 批 复 了 《岩溶 地 区 石 漠 化 综合 治理 规 
划 大 纲 (2006 一 2015) 》, 在 西南 8 个 省 (区 直辖 市 ) 100 个 石 漠 化 严重 县 启动 看 漠 化 综合 治理 工程 试点 ,探索 
石 漠 化 综合 治理 的 模式 与 途径 。2014 年 已 扩大 至 316 个 县 , 占 到 全 国 455M 锥 石 漠 化 县 的 69.5% 。 截 止 2015 
年 底 , 石 漠 化 治理 重点 工程 县 投入 中 央 预 算 内 专项 资金 119 亿 元 , 石 漠 化 在 地 退化 趋势 得 到 有 效 遏 制 , 石 漠 化 
土地 面积 已 实现 由 持续 增加 向 “ 净 减 少 ” 的 重大 转变 。 然 而 ,由 对 咯 斯 特地 质 背 景 的 复杂 性 及 石 漠 化 生态 恢 
复 过 程 中 的 不 确定 性 ,喀斯特 地 区 石 漠 化 治理 成 效 在 空间 H 具 有 明显 差异 ''"|, 据 2012 年 6 月 国家 林业 局 公 
布 的 《中 国 石 漠 化 状况 公报 》, 与 2005 年 相 比 ,广西 石 漠 先 土 地 面积 减少 最 多 ,贵州 次 之 ,云南 最 少 。 石 汉化 
治理 的 核心 是 植被 恢复 ,因此 ,研究 西南 喀斯特 地 区 县 域 植 被 恢复 变化 ,对 开展 喀斯特 石 漠 化 治理 工程 综合 成 
效 评估 ,改进 喀斯特 地 区 石 漠 化 治理 措施 与 后 续 玉 程 规划 具有 重要 现实 意义 。 

现 有 研究 已 揭示 了 喀斯特 地 区 植被 增加 的 医 域 差异 比较 明显 ,在 年 际 尺 度 上 ,西南 喀斯特 地 区 植被 变化 
与 温度 和 降水 均 无 显著 相关 "” ,在 区 域 尺度 了 总 气候 因子 (年 均 降水 和 年 均 温 度 ) 不 是 近 12 年 来 喀斯特 生态 
恢复 好 转 的 主导 因素 。 生 态 工 程 的 实施 是 导致 喀斯特 石 漠 化 地 区 生态 恢复 好 转 的 原因 。 然 而 ,大 区 域 
尺度 ,特别 是 作为 我 国 石 漠 化 集中 分 布 的 汗 柱 黔 三 省 ,喀斯特 地 貌 类 型 复杂 多 样 ,生态 建设 工程 实施 以 来 植被 
恢复 状况 如 何 ,不 同 喀斯特 地 貌 类 型 区 植被 变化 的 主导 因素 存在 怎样 的 区 域 差 异 目前 尚未 明确 ,因此 , 急 
需 开 展 西 南 喀斯特 地 网 植被 变化 主导 因素 空间 非 平稳 性 研究 。 

根据 “地 理学 第 -定律 ”'5 ,普通 的 全 局 模型 使 用 的 假设 前 提 是 数据 本 身 在 统计 上 是 独立 的 ,并 呈正 态 
分 布 ” ,忽略 了 研究 问题 的 空间 效应 , 且 估 计 的 回归 系数 是 一 个 常数 ,将 较 难 识别 区 域 尺度 上 因 地 理 位 置 的 
变化 而 引起 的 喀斯特 植被 变化 及 其 影响 因子 间 关 系 的 空间 非 平 稳 性 。 考 虑 样本 空间 地 理 位置 的 GWR 模型 
有 效 改 进 了 宝 量 闻 相 互 关系 的 空间 非 平稳 性 分 析 , 能 有 效 克服 喀斯特 地 区 地 理 要 素 具有 显著 的 空间 异 质 性 和 
空间 关联 性 的 难题 。 因 此 ,本 研究 基于 GWR 模型 ,利用 长 时 间 序 列 遥 感 数据 分 析 近 30 年 来 (1982 一 2011 
年 ) , 汗 桂 黔 三 省 植被 变化 趋势 及 其 主导 因素 ,解析 县 域 尺度 植被 变化 的 空间 效应 与 集聚 性 规律 ,揭示 植被 变 
化 主导 因素 的 空间 非 平稳 性 与 区 域 差 异 ,以 期 为 后 续 石 漠 化 治理 工程 的 实施 与 工程 成 效 评估 提供 科学 依据 。 


1 研究 地 区 与 研究 方法 


1.1 研究 区 概况 

本 研究 区 域 位 于 我 国 西南 喀斯特 地 区 的 溃 桂 上 黔 三 省 (97.81° 一 111.67°E,21.53° 一 28.93°N) ,是 我 国 喀 斯 
特地 貌 发 育 的 集中 分 布 区 "" 。 该 区 域 气候 温和 湿润 ,年 均 降水 量 1000 一 1600mm ,降雨 时 空 分 布 不 均 , 干 湿 交 
百 频 繁 (1 中, 属 亚热带 季风 气候 区 。 溃 桂 黑 地 区 地 带 性 植被 为 中 亚热带 常 绿 阔 叶 林 ,植被 类 型 多 样 ,主要 类 
型 有 常 绿 针 叶 林 、 常 绿 闭 叶 林 、 灌 从 及 灌 草 从 等 。 脆 弱 的 地 质 背 景 加 上 长 期 以 来 人 类 的 不 合理 干扰 ,使 研究 区 
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内 自然 植被 遭 到 破坏 ,土壤 侵蚀 严重 , 基 岩 大 面积 裸露 , 石 漠 化 问题 突出 。 截 止 2011 年 底 , 汗 桂 黔 三 省 仍 有 石 
漠 化 土地 面积 779.0 万 hm , 占 我 国 石 漠 化 土地 总 面积 的 64.9% , 主要 分 布 于 峰 林 平原 . 峰 从 洼地、 岩溶 高 原 、 
岩溶 峡谷 . 断 陷 盆 地 等 喀斯特 地 瑶 类 型 区 (图 1)。 
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图 1 研究 区 概况 


Fig.l General Situation of the research area 


1.2 ”数据 来 源 与 处 理 
1.2.1 NDVI 数据 的 获取 与 处 理 

本 文 NDVI 数据 采用 的 是 CIMMS NDVI38g 数据 集 产品 , 空 间 分 辨 率 是 8km , 时 间 分 辨 率 15d, 从 1982 一 
2011 年 共 792 期 影像 。 该 数据 集 对 以 前 的 版 本 进行 了 改进 与 延伸 ,运用 改进 的 云 掩 膜 技 术 对 与 植被 绿 度 无 
关 的 因子 进行 了 校正 呈 ; 。 为 减少 秀和 大 气 造 成 的 噪音 ,将 一 月 两 次 数据 中 的 最 大 值 作为 当月 NDVI, 并 考虑 
喀斯特 地 区 植被 物候 特征 ,将 每 年 4 一 正月 的 平均 值 作 为 当年 GSN。 以 2001 年 作为 生态 工程 实施 节点 六 ， 
1982 一 2000 年 为 工程 实施 之 前 植被 状况 ,2001 一 2011 年 为 工程 实施 后 的 植被 状况 ,两 个 时 间 段 的 差 值 用 来 表 
征 喀斯特 地 区 植被 变化 情况 ,并进 行 归 一 化 处 理 。 
1.2.2 植被 变化 影响 因素 指标 体系 获取 与 处 理 

为 了 保证 选取 的 影响 因素 的 全 面 性 以 及 县 域 尺 度 上 数据 的 可 获取 性 ,本 研究 选择 了 与 植被 生长 关系 密切 
的 自然 因素 (万 员 热 、 区. 地质 背 景 等 ) 和 人 类 活动 因素 等 , 共 17 个 指标 。 具 体 数 据 来 源 如 下 :气象 数据 来 自 中 
国 气象 数据 网 ( http://data.cma.cn/) ,包括 汗 桂 黔 地 区 的 2001 一 2011 年 逐年 年 均 气 温 、 降 水 .湿度 指数 .干燥 
度 .积温 等 本 研究 用 县 域 尺度 上 11a 的 气象 因子 均值 表征 生态 工程 实施 后 该 县 的 平均 气候 条 件 。2005 年 
和 2010i 年 八 口 密度 .CDP .土地 利用 类 型 坡度 不 同 植被 和 土壤 条 件 下 土壤 有 机 碳 等 数据 从 中 国 科 学 院 资源 
环境 科学 数据 中 心 (http://www.resdc.cn/) 下 载 ; 土 层 厚 度数 据 来 源 于 寒 区 旱 区 科学 数据 中 心 (http://westdc. 
westgis.ac.cn) ;退耕 还 林 工 程 资 金 投入 和 造林 面积 数据 由 汗 桂 黔 3 省 林业 厅 提 供 。 退 耕 还 林 工 程 投入 包括 退 
耕 还 林 投入 .荒山 造林 投入 和 封山育林 投入 及 现金 补助 4 个 部 分 。 本 研究 中 生态 工程 投入 和 造林 面积 为 各 县 
从 2001 一 2011 年 的 单位 面积 上 的 资金 总 投入 及 造林 总 面积 。 为 确保 计量 单位 的 统一 性 及 各 因子 间 的 可 比 
性 ,本 研究 对 以 上 各 因素 进行 了 归 一 化 处 理 。 
1.3 植被 变化 主导 因素 识别 方法 
1.3.1 空间 自 相 关 性 分 析 

变量 存在 空间 自 相 关 性 是 应 用 GWR 模型 分 析 的 前 提 条 件 。GSN 及 其 影响 因素 属于 地 理 空 间 数 据 , 具 有 
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空间 自 相关 性 和 空间 非 平稳 性 等 特征 ,因而 它们 的 空间 关系 可 用 全 局 空间 自 相关 莫 兰 指数 ( 1 ) 和 局 部 空间 
自 相 关 莫 兰 指数 ( 厂 ) 进行 表征 。 
1.3.2 多重 共 线 性 分 析 

进行 GWR 模型 回归 时 ,如 果 变 量 间 相关 性 程度 高 ,出 现 信息 重 和 至 现象 ,就 容易 产生 多 重 共 线 性 问题 。 相 
关 分 析 是 判断 多 重 共 线性 问题 .剔除 元 余 变量 的 基本 方法 。 如 果 相 关系 数值 超过 0.8, 则 变量 之 间 存 在 较 强 的 
相关 关系 :5 ,相关 系数 小 于 0.8, 则 所 选 因子 之 间 不 存在 多 重 共 线性 问题 。 本 研究 利用 SPSS 软件 中 主 成 分 分 
析 方 法 (PCA) 对 原始 影响 因子 数据 进行 变换 和 处 理 , 使 原来 多 个 相互 影响 的 变量 转换 成 几 个 互 不 相 产 的 综 
合 变量 , 既 保 留 了 原始 影响 因素 的 主要 信息 ,又 消除 了 变量 间 的 多 重 共 线 性 ,得 到 影响 植被 变化 的 综合 主导 因 
素 !5 。 本 研究 用 方差 最 大 旋转 矩阵 将 主 成 分 中 各 变量 的 贡献 率 向 最 大 和 最 小 转变 ,使 每 个 主 成 分 上 具有 最 
高 载荷 的 因子 数 最 少 ,从 而 使 得 对 主 成 分 的 解释 更 清晰 5 。 
1.3.3 ”地理 加 权 回 归 模 型 

空间 平稳 是 指 与 地 理 空间 关联 的 自 变 量 与 因 变 量 之 间 不 存在 空间 差异 , 即 B,(wh ao) 不 随 自 变量 x, 位 置 
而 变化 ,反之 则 为 空间 不 平稳 '” 。 空 间 数据 一 般 具 有 空间 非 平稳 性 的 特征 ,GWR 模型 侈 许 在 不 同 的 地 理 空 
间 有 不 同 的 空间 关系 存在 ,因此 它 是 指 一 种 相对 简单 又 有 效 的 探测 空间 非 平稳 性 的 新 方法 局] 。GWR 模型 最 
早 是 由 Brunsdon 等 人 提出 中 ,用 于 研究 自 变量 在 空间 上 的 局 部 非 平稳 性 其 具体 公式 如 下 : 

y; = > Bi(uiyv) xi Te， 

式 中 , Bi 指 局 部 系数 , (u,v;) 为 地 理 位 置 , x, ; 为 局 部 因素 , 即 6 个 主要 因素 , se, 指 的 是 服从 均值 为 零 的 独立 
正 态 分 布 的 误差 。 对 空间 权重 矩阵 进行 计算 的 权重 函数 选择 的 是 高 斯 模型 , 核 函 数 是 自 适 应 性 双 平 方 方 法 ， 
最 优 带宽 选择 利用 的 高 斯 方法 ,利用 AICc 信息 准则 法 确定 最 做人 带宽 。 将 2001 一 2011 年 4 一 11 月 平均 GSN 
与 1982 一 2000 年 平均 GSN 的 差 值 作为 GWR 模型 的 因 变 量 4 利 用 PCA 提取 主 成 分 作为 自 变 量 ,得 到 每 个 县 
域 的 影响 因素 的 回归 系数 ,分 析 单 因素 的 空间 非 平 稳 性 。 回 归 系 数 绝对 值 越 大 ,该 因素 对 植被 变化 的 影响 越 
大 ;反之 , 越 小 ;回归 系数 等 于 0, 两 者 没有 关系 个 同时 利用 最 小 二 乘 方 法 进行 全 局 回归 分 析 , 并 对 比 两 模型 的 
结果 。 为 了 进一步 验证 各 回归 系数 在 空间 上 大 浊 具有 非 平稳 性 ,对 各 主 成 分 的 回归 系数 进行 空间 自 相关 分 
析 。 将 每 个 县 域 回 归 系 数 绝对 值 最 太 和 值 对 应 的 影响 因素 作为 影响 本 县 植被 恢复 的 主导 因素 ,进一步 分 析 主 导 
因素 的 空间 差异 性 。 


2 结果 分 析 


2.1 生态 工程 前 后 西南 所 省 植被 覆盖 变化 总 体 分 析 

根据 西南 喀斯特 地 区 重大 生态 工程 实施 背景 及 Tong 等 中 对 工程 实施 后 植被 履 盖 变化 的 突变 检验 分 析 ， 
将 西南 三 省 植被 覆盖 时 间 变 化 分 为 工程 实施 前 (1982 一 2000 年 ) 与 实施 后 (2001 一 2011 年 )。 相 较 于 工程 实 
施 前 , 汗 桂 黔 县 域 尺度 平均 GSN 在 工程 实施 后 均 有 所 增加 ,但 增加 程度 在 空间 上 具有 空间 差异 (图 2) 。 空 间 
变化 上 在 省 域 尽 度 , 工 程 实施 前 后 汗 桂 黔 GSN 平均 增加 了 0.035 ,其 中 贵州 和 广西 GSN 分 别 相对 增加 了 
0.045 和 ,0.045 ,云南 GSN 相对 增加 较 少 (0.024) ;生态 工程 前 后 溃 桂 黔 三 省 各 县 平均 GSN 变化 值 范 围 介 于 
0.0002 和 0.07 之 间 , 且 变异 系数 为 43% ,说 明 各 县 实施 生态 工程 后 都 取得 成 效 , 但 各 县 成 效 有 差异 。 

对 比 不 同 石 漠 化 治理 工程 类 型 区 GSN 变化 发 现 ,不 同 石 漠 化 治理 工程 类 型 区 GSN 呈 总 体 增 加 趋势 (图 
3) ,植被 增加 的 区 域 占 整个 研究 区 域 的 54.3% 。 不 同 工 程 类 型 区 植被 变化 存在 一 定 的 差异 :岩溶 区 植被 增加 
面积 占 岩 游 区 总 面积 的 55.49% ;工程 区 中 岩溶 村 谷 . 峰 林 平原 岩溶 高 原 . 峰 从 洼地 等 以 植被 增加 为 主 , 分 别 
占 各 工程 区 总 面积 的 83.00% ,78.14% ,69.64% ,56.76% , 其 中 ,岩溶 覃 和 谷 植被 增加 比例 最 多 ;中 高 山 、 岩溶 峡 
谷 .岩溶 断 陷 盆地 以 植被 不 变 为 主 。 上 述 结 果 说 明 近 30 年 来 滇 桂 黔 3 省 植被 总 体 呈 增加 趋势 ,生态 工程 的 实 
施加 快 了 西南 三 省 植被 覆盖 的 恢复 速率 ,但 不 同 区域 特别 是 不 同 喀斯特 地 貌 类 型 区 植被 覆盖 变化 存在 显著 
的 空间 差异 性 。 
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图 2 生态 工程 实施 前 后 NDVI 变化 情况 
Fig.2 Changes of NDVI before and after the implementation of ecological engineering projects 


2.2 植被 覆盖 变化 的 空间 格局 差异 性 分 析 1.6 


GSN 空间 自 相 关 性 检验 分 析 显 示人 研究 区 7 的 正 态 人 记 4 
统计 量 4 值 为 30.30, 大 于 正 态 分 布 函数 在 0.01 水 平 RIG_ 这 1 
的 临界 值 2.58 , 7 为 0.90 ,表明 研究 区 GSN 具有 最 著 的 十 08 
正 空间 自 相 关 性 。 从 7, 空间 分 布 ( 图 4) 可 以 看 旧部 分 E 
岩溶 模 谷 . 峰 从 洼地 和 峰 林 平原 地 区 呈现 植被 寿 数 高 值 三，， 
集聚 ,是 植被 密集 分 布 区 ; 央 溶 峡谷 、. 兰 浴 断 陷 盆 地 区 的 Ce 
植被 指数 则 呈现 低 值 集聚 , 即 植 袖 履 疹 较 小 。 因 此 , 工 六 年 问 可 息 答 久 
程 实施 后 汗 桂 黔 3 省 平均 GSN 在 整体 上 具有 显著 的 空 稚 EE 
扫 


间 差 异性 ,在 局 部 区 域 上 具有 显著 的 空间 集聚 性 。 
2.3 演 桂 黔 植被 变化 的 主要 因素 分 析 

吻 除 与 植被 指数 赠 长 值 不 显著 相关 的 变量 ,利用 0 
PCA 方法 处 理 剩余 自 克 量 之 问 的 多 重 共 线 性 ,得 到 影 Fig.3 The change of NDVI in different ecological engineering 
响 因 素 的 主 威 分 分 析 结 果 和 各 成 分 的 得 分 系数 矩阵 , 利 
用 旋转 成 分 矩阵 ( 表 1) 及 相关 分 析 , 可 以 更 清晰 的 查看 每 个 因素 对 各 主 成 分 的 贡献 率 。 每 个 主 成 分 可 以 根据 
载荷 较 关 的 影响 因素 将 主 成 分 进行 分 类 ( 表 2) 。 在 本 研究 结果 显示 ,前 6 个 主 成 分 (PCA1 一 PCA6) 的 累计 贡 
献 率 为 78.33% ,可 较 好 地 表征 原 数 据 '” 。 以 植被 指数 增长 值 作为 因 变量 ,与 6 个 主 成 分 进行 OLS 回归 ( 表 
2) 结果 表 明 6 个 主 成 分 的 方差 膨胀 因子 都 为 1.00, 则 6 主 成 分 之 间 不 存在 多 重 共 线 性 。 但 对 OLS 多 元 回归 
的 残 差 进行 空间 自 相关 性 检验 , 7, 为 0.53 ,2 统计 值 为 17.90, 说 明 残 差 在 空间 上 都 存在 显著 地 正 空间 自 相关 
性 ,因此 需要 利用 GWR 模型 进一步 分 析 影 响 因素 的 空间 差异 性 。 
2.4 不 同 影响 因素 对 植被 覆盖 变化 的 空间 非 平稳 性 
2.4.1 植被 变化 主要 影响 因素 的 GWR 分 析 

以 GSN 增长 值 为 因 变 量 ,以 对 6 个 主 成 分 作为 自 变 量 进行 空间 计量 GWR 分 析 , 发 现 残 差 莫 兰 指数 与 
OLS 模型 相 比 有 所 降低 (元 为 0.11,2(7) 为 3.66)。 同 时 结合 相关 性 和 矩阵 和 主 成 分 定性 分 析 ( 表 3) ,结果 显示 


生态 工程 区 Ecological area 


area 
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4 2001 一 2011 年 平均 生长 季 NDVI 增长 值 的 局 部 空 


因 琵 ,利用 GWR 模型 对 喀斯特 地 区 植被 变 


Fig.4 The local spatial autocorrelation of added average gr 


CWR 模型 解释 了 GSN 变化 的 79% , 远 优 于 全 局 模型 (R?=0083)，。 


化 的 影响 因素 进行 回归 分 析 具 有 合理 性 


oO 


表 1 旋 


Table 1 ation component matrix 


O] 


Bl- | 


VI in 2001 一 2011 


4 


各 因子 贡献 率 6 个 主 成 分 Six principal components 

Contribution rate of each factor 2 ie 4 5 6 
资金 投入 Capital input -0.02 0.12 0.88 0.08 0.14 0.04 
造林 面积 Afforestation area -0.24 0.18 0.87 -0.07 —0.08 —0.08 
不 同 土壤 条 件 下 SOC SOC on soil type -0.42 0.19 —0.03 0.29 0.60 -0.17 
不 同 植被 条 件 下 SOC SOC on grass f5 0.32 -0.06 -0.02 -0.07 0.74 0.18 
土壤 氮 储 量 Soil nitrogen stora —0.64 -0.06 -0.03 0.40 0.17 0.17 
>0% 积温 >0%C accumulated temperature 0.87 -0.27 -0.13 0.13 0.05 0.24 
气温 Temperature ” 轧 0.85 -0.25 -0.29 0.06 -0.03 0.10 
>10C 积 温 >103GEEecumulated temperature 0.88 -0.27 -0.15 0.12 0.04 0.22 
Wor hsity 0.04 -0.29 0.47 -0.04 -0.01 0.57 
0.14 0.25 -0.16 0.09 0.03 0.73 
-0.26 0.51 -0.07 0.08 0.02 -0.67 

0.27 —0.85 -0.21 0.09 -0.20 0.020 

粒 含 量 Sand content -0.22 0.90 0.08 0.00 -0.21 0.00 
粉 粒 含量 Silt content -0.10 -0.06 0.25 -0.17 0.79 —0.04 
干燥 度 Aridity -0.19 0.18 -0.28 0.06 -0.03 0.10 
降水 量 Precipitation 0.65 -0.11 -0.25 0.58 0.10 0.22 
湿度 指数 Humidity index —0.08 0.11 -0.01 0.89 -0.11 —0.04 


提取 方法 : 主 成 份 ; 旋转 法 :具有 Kaiser 标准 化 的 正 交 旋 转 法 ; a. 旋转 在 6 次 迭代 后 收敛 
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表 2 基于 主 成 分 的 OLS 模型 参数 估计 
Table 2 PCA based Estimation Parameter Values of OLS 


变量 参数 估计 值 标准 误差 , 方差 膨胀 因子 主导 因子 
Variables Parameter estimation Standard error VIF[c] Dominant factors 
常数 项 Constant term 0.51 0.01 44.44 基础 自然 条 件 
PCAI1 0.06 0.01 5.52 1.00 气候 因子 
PCA 2 0.00 0.01 —0.05 1.00 土壤 质地 
PCA 3 0.02 0.01 1.85 1.00 人 类 活动 
PCA 4 0.03 0.01 2.30 1.00 水 分 有 效 性 
PCA 5 0.07 0.01 6.39 1.00 土壤 养分 
PCA 6 0.03 0.01 2.82 1.00 社会 经 济 移 件 


表 3 GWR 模型 与 全 局 回归 模型 拟 合 优 度 


Table 3 ” Comparison of coefficients of determination for the GWR and Least squared regression models 


模型 Model 拟 合 优 度 R? 调整 拟 合 优 度 Adjusted R? 局 部 拟 合 优 度 Local R? 
地 理 加 权 回 归 模 型 CWR 0.79 0.71 0.25 一 0.81 
全 局 回归 模型 OLS 0.23 0.22 


为 了 进一步 验证 各 主 成 分 回归 系数 在 空间 上 是 否 平 稳 ,对 回归 系数 进行 空间 自 相关 性 分 析 ( 表 4) ,得 出 
GWR 模型 的 各 回归 系数 的 二 均 大 于 0,2Z(1) 均 大 于 2.58 ,说 明 在 此 % 显 著 水 平 下 ,各 项 系数 具有 显著 地 正 空 间 
自 相 关 性 , 即 这 些 回 归 系 数 在 空间 上 具有 空间 非 平稳 性 。 因 此 ,利用 CCWR 模型 可 以 较 详细 得 考虑 不 同 影 响 
因素 在 不 同 的 地 理 位 置 对 植被 变化 的 定量 影响 。 


表 4 各 因素 回归 系数 的 莫 兰 指数 及 Z 统计 量 


Table 4 Moran's I and Z,Scores of regression coefficients for each factor 


影响 因子 第 一 主 成 分 第 二 主 成 分 第 三 主 成 分 第 四 主 成 分 第 五 主 成 分 第 六 主 成 分 
Impacting factors PCA1 PCA2 PCA3 PCA4 PCAS PCA6 
Moran's J 0.68 0251 0.42 0.57 0.76 0.41 
区 23.04 17.21 14.68 19.10 25.50 14.31 


2.4.2 ” 单 因素 对 植被 变化 影响 的 罕 间 非 平稳 性 分 析 

研究 区 内 各 主 成 分 的 回归 又 数 空间 分 布 差异 较 大 ( 图 5) ,各 成 分 的 回归 系数 既 有 正 值 ,又 有 负 值 , 且 占 样 
本 总 数 的 比例 都 较 高 ,表明 该 区 域 各 影响 因素 与 植被 变化 的 关系 较 复杂 ,各 影响 因子 对 植被 变化 影响 程度 具 
有 显著 的 空间 异 质 性 。 气候 因素 对 植被 作用 程度 的 空间 分 布 具 有 一 致 性 ,大 致 规律 是 沿 云南 高 原 、 黔 西山 地 
高 原 . 黔 北山 呈 仿 高原 斜 坡地 一 广西 贫 地 的 海拔 梯度 呈现 递减 特征 ,同时 向 东 至 峰 林 平原 地 区 递增 ,体现 了 喀 
斯 特地 区 水 热 空间 分 配 差 异 明 显 ,致使 植被 对 水 热 的 响应 也 具有 空间 差异 。 土 壤 质 地 对 植被 的 影响 有 两 个 分 
布地 带 :在 过 贵 鹿 原 地 区 大 致 沿 长 江 流域 和 珠江 流域 的 分 水 岭 向 两 边 递减 , 呈 对 称 分 布 ;在 两 广 丘陵 地 区 大 致 
游 启 西 岳 陵 向 广西 盆地 递减 ,体现 了 河流 上 下 游 的 土壤 质地 差异 。 第 五 主 成 分 全 加 了 土壤 养分 的 影响 ,在 涯 
次 区 域 竺 壤 属 性 与 植被 增长 以 正 相 关 关 系 为 主 , 知 深 峡谷 和 峰 林 平 原 东 部 区 域 以 负 相 关 关 系 为 主 。 

人 类 活动 对 植被 具有 双重 作用 ,表现 在 为 恢复 生态 环境 实施 的 一 系列 生态 保护 与 建设 工程 措施 的 正 向 作 
用 和 为 发 展 社会 经 济 水 平 进行 的 一 系列 不 合理 的 生产 活动 的 负面 影响 。 在 断 陷 倪 地 、 峰 从 洼地 西部 及 峰 林 平 
原 东 部 生态 工程 对 植被 恢复 起 到 了 促进 作用 , 知 溶 模 谷 峡谷、 高 原 和 部 分 峰 从 洼地 、 峰 林 平原 区 人 类 活动 的 
逆向 干扰 抑制 了 植被 增长 。 社 会 经 济 耦合 地 形 因素 对 植被 生长 的 影响 空间 复杂 性 更 高 ,大 致 沿 云贵 高 原 向 广 
西 盆地 过 渡 的 斜坡 地 带 与 植被 增长 呈正 相关 关系 ,并 向 东南 方向 的 广西 盆地 区 域 过 渡 为 负 相关 关系 。 这 反映 
了 人 类 在 地 形 复杂 的 斜坡 地 带 进行 的 不 合理 的 生产 活动 对 植被 增长 产生 逆向 干扰 。 
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5 气候 因子 (a) 土壤 质地 (b) 、 人 类 活动 (c) 、 水 分 有 效 性 (d) .土壤 养分 (e) 和 社会 经 济 条 件 (f) 等 因素 对 GSN 变化 的 空间 非 平 稳 性 


Fig.5 The spatial non-stationary of climate( a) soil texture (b ) human activities ( c ) moisture availability ( d) ,soil nutrients ( e ) social 


economical condition(f) on GSN at county level 
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图 6 主导 因素 对 植被 变化 影响 的 区 域 差异 


Fig.6 Regional difference of dominant factors on vegetation dynamics 


2.4.3 植被 变化 主导 因素 的 空间 格局 分 析 

为 了 进一步 分 析 生态 工程 背景 下 具有 复杂 地 质 背景 的 路 斯 特地 区 造成 植被 恢复 差异 的 原因 ,本 研究 综合 
以 上 分 析 的 单个 影响 因素 的 空间 非 平 稳 性 特征 ,利用 MATEAB 编程 提取 每 个 县 域 对 植被 生长 影响 最 大 的 因 
素 (图 6) ,发 现在 生态 工程 背景 下 影响 植被 变化 的 主导 因素 随 不 同 喀斯特 地 形 地 貌 特 征 的 变化 而 变化 , 且 各 
影响 因子 与 植被 变化 既 有 正 相关 关系 又 有 负 相 关 闫 系 。 人 类 活动 和 社会 经 济 条 件 是 峰 林 平原 区 . 断 陷 盆 地 区 
及 汗 桂 黔 交界 的 石 漠 化 严重 地 区 的 植被 变化 的 训导 因素 ,而 气候 因素 .土壤 养分 ,水 分 有 效 性 等 自然 因素 对 岩 
洲 峡 谷 区 , 槽 谷 区 .高原 区 及 部 分 峰 傣 洼地 区 的 植被 变化 影响 均 较 为 显著 。 其 中 , 峰 林 平原 区 植被 变化 受 社会 
经 济 条 件 的 影响 较 大 , 峰 丛 洼地 区 党 气候 因素 及 人 类 活动 的 共同 作用 (部 分 地 区 植被 变化 也 受 水 分 有 效 性 及 
土壤 养分 的 影响 ) ,喀斯特 高 厌 区 植被 变化 受气 候 因子 .土壤 养分 和 水 分 有 效 性 等 自然 因素 的 影响 较 大 ,中 高 
山地 区 的 植被 主要 受降 水 及 状 类 活动 的 影响 ;而 岩溶 峡谷 区 植被 变化 的 影响 以 气候 因子 为 主导 , 断 陷 倪 地 区 
植被 主要 受 人 类 活动 及 社会 经 济 条 件 的 影响 ,岩溶 模 谷 区 植被 变化 受 水 分 有 效 性 影响 显著 。 因 此 ,影响 西南 
喀斯特 地 区 植被 变化 的 主导 因素 在 空间 上 具有 显著 的 区 域 差 异性 ,不 同 喀斯特 地 貌 类 型 间 植 被 变化 主导 因素 
差异 较 大 。 


西南 喀斯特 区 域 具有 复杂 的 地 质 背 景 ,高 度 的 景观 异 质 性 ,复杂 的 二 元 水 文系 统 , 特 有 的 植被 种 类 和 高 强 
度 的 类 类 活动 ,这些 复 杂 的 环境 背景 决定 了 植被 变化 的 空间 差异 性 和 影响 植被 变化 的 影响 因子 的 多 样 性 和 复 
杂 性 4 本 研究 揭示 了 受 地质 背 景 制约 的 西南 喀斯特 地 区 植被 在 气候 变化 和 人 类 活动 等 作用 下 植被 变化 的 主 
导 因 素 及 由 于 区 域 背 景 条 件 差异 导致 的 主导 因素 变化 的 空间 非 平稳 性 。 

从 全 局 看 , 黔 东 桂 东南 及 云南 地 区 的 岩溶 断 陷 贫 地 等 地 区 ,主要 受气 候 因 素 的 影响 ,热量 因子 与 植被 覆 
盖 呈 负 相 关 关 系 。 近 年 来 喀斯特 地 区 植被 覆盖 变化 呈 上 升 趋势 ,但 2001 一 2011 年 西南 喀斯特 地 区 呈 干 旱 化 
趋势 ,而 该 区 年 降水 量 在 800mm 以 上 ,所 以 降水 减少 并 不 能 导致 植被 发 生 干 旱 胁迫 ,喀斯特 地 区 干旱 化 
可 能 是 由 于 气温 升 高 导致 蒸 散 增加 所 致 ,这 与 王 兆 礼 等 .研究 的 中 国 西南 喀斯特 区 域 NDVI 与 SPEI 呈 负 相 
关 结论 一 致 。 云 贵 高 原 东 南部 及 其 与 广西 盆地 的 过 渡 地 带 石 漠 化 问题 严重 ,地 表 植 被 遭 到 严重 破坏 , 随 着 大 
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规模 生态 工程 的 实施 ,植被 覆盖 逐渐 增加 。 但 由 于 实施 时 间 相 对 较 短 ,大 部 分 恢复 区 域 多 处 于 植被 恢复 初期 ， 
以 草 灌 为 主 ,加 上 复杂 的 地 质 背 景 ,植被 覆盖 变化 影响 因素 复杂 ,受气 候 因 素 .社会 经 济 条 件 .土壤 状况 及 人 类 
活动 干扰 的 共同 作用 (图 6) 。 

不 同 喀斯特 地 貌 类 型 条 件 下 ,植被 变化 的 主导 因素 存在 显著 的 空间 非 平稳 性 ,与 不 同 地 貌 类 型 的 水 文 地 
质 特征 密切 相关 。 峰 林 平 原 区 (图 61 ) 地 形 相对 较为 平坦 交通 条 件 好 、 社 会 经 济 活动 频繁 ,植被 变化 受 社会 
经 济 条 件 的 影响 较 大 。 峰 从 洼地 区 (图 6 了) 水 热 条 件 相 对 较 好 ,生态 工程 实施 以 来 植被 增长 显著 ,但 仍 存在 
不 合理 人 类 活动 干扰 ,植被 变化 受气 候 因 素 及 人 类 活动 的 共同 影响 显著 ;同时 ,由 于 该 区 土地 资源 分 散 和 典型 
二 元 水 文 地 质 结 构 ,地 表 水 系 缺 乏 ,部 分 地 区 植被 变化 也 受 水 分 有 效 性 及 土壤 养分 的 影响 。 喀 斯 特 高 原 区 
(图 6VI) 云雾 多 日照 少 ,太阳 辐射 能 量 低 ,海拔 落差 较 大 ,该 区 植被 变化 受气 候 因 子 .土壤 养分 和 水 分 有 效 性 
等 自然 因素 的 影响 较 大 。 中 高 山地 区 (图 6JV ) 山高 坡 陡 ,自然 条 件 较 差 ,局 部 水 资源 匮乏 ,植被 变化 主要 受降 
水 及 人 类 活动 的 影响 。 岩 溶 峡 谷 区 (图 6V ) 以 典型 的 中 亚热带 山区 气候 和 南亚 热带 干 热 河谷 为 特征 ,植被 变 
化 的 影响 以 气候 因子 为 主导 。 断 陷 盆 地 区 的 盆地 内 (图 6 碍 ) 地 形 平坦 、 光 热 条 件 好 , 典 类 活动 频繁 , 而 四 周 为 
山地 、 地 形变 化 大 ,植被 变化 主要 受 人 类 活动 及 社会 经 济 条 件 的 影响 。 岩 溶 模 舍 区 人 图 6 班 ) 具有 一 定 厚 度 的 
土 层 分 布 ,水 土 资源 条 件 较 好 ,土地 承载 能 力 相对 较 高 ,但 该 区 地 质 构 造 复杂 ,第 出 现 * 苹 挂 式 ” 地 下 河 ,地 表 
水 系 缺乏 ,植被 变化 受 水 分 有 效 性 影响 显著 。 

人 类 活动 对 植被 增长 具有 双重 作用 ,并 通过 土地 利用 方式 反映 对 生 容 环境 的 影响 程度 。 一 方面 ,喀斯特 
地 区 石 漠 化 是 在 典型 的 生态 环境 脆弱 区 ,高 强度 人 类 活动 干扰 导致 喀斯特 植被 破坏 、 其 至 土地 石 漠 化 的 现 
象  。 另 一 方面 ,生态 保护 与 建设 工程 措施 促进 了 植被 覆盖 的 增加 ,Tong 等 六 剔除 气象 因素 对 植被 变化 的 
影响 分 析 了 黔 滨 桂 地 区 人 类 活动 对 植被 的 干扰 , 发现 生态 工程 对 植被 恢复 的 正 向 作用 。 然 而 ,由 于 不 同 地 理 
位 置 自 然 条 件 、. 地 质 背 景 和 干扰 方式 及 强度 等 方面 的 差异 ,退耕 还 林 、 石 漠 化 治理 等 生态 工程 成 效 具有 空间 差 


异性 , 李 吴 等 ”对 贵州 毕节 地 区 的 退耕 还 林 工 程 效果 进行 评 佑 发 现 退 耕 还 林 等 生态 工程 的 实施 加 快 了 喀 斯 
特地 区 植被 恢复 和 生态 重建 ,但 赫 章 ,大方 等 喀斯特 地 区 在 巨大 的 人 口 压力 下 , 仍 存在 较为 严重 的 植被 退化 现 
象 。 因 此 ,喀斯特 地 区 生态 建设 与 保护 是 项 长 期 而 艰巨 的 任务 ,综合 考虑 植被 变化 主导 因素 的 区 域 差 异 来 
调整 或 改进 不 同 喀斯特 区 域 背景 条 件 下 失态 工程 的 措施 ,促进 喀斯特 地 区 生态 环境 的 恢复 。 


4 结论 


本 研究 采用 长 时 间 序 列 遂 感 数据 ,综合 运用 空间 自 相 关 分 析 、PCA 和 GWR 等 研究 方法 ,分 析 生 态 工程 实 
施 以 来 汗 桂 黑 喀 斯 特 植被 变化 及 其 主要 影响 因素 的 空间 非 平 稳 性 。 研 究 表 明 近 30 年 来 演 桂 黔 三 省 植被 总 体 
呈 增 加 趋势 ,生态 工程 的 实施 加 快 了 西南 三 省 植被 覆盖 的 恢复 速率 ,但 不 同 区 域 .特别 是 不 同 喀斯特 地 貌 类 型 
区 植被 覆盖 变化 存在 器 著 的 空间 差异 性 。 植 被 变化 的 空间 差异 性 主要 是 气候 土壤、 地形 和 人 类 活动 等 因素 
引起 的 ,影响 植被 变 借 的 主导 因素 随地 理 位置 的 变化 而 不 同 , 具 有 显著 的 空间 非 平 稳 性 。 

喀斯特 地 区 自然 条 件 优越 ,但 复杂 的 坡度 等 地 形 条 件 炙 加 高 强度 的 人 类 活动 干扰 使 生态 环境 严重 退化 。 
由 于 其 独特 的 地 质 背景 ,生态 工程 实施 难度 大 。 因 此 ,加 大 生态 工程 投入 ,完善 生态 工程 措施 ,在 具体 实施 过 
程 电 要 综合 考虑 影响 植 侦 生 长 的 影响 因素 以 及 不 同 地 貌 类 型 区 主导 因素 的 空间 差异 ,对 不 同 工 程 区 进行 科学 
规划 与 优化 分 区 ,因地制宜 采取 生态 移民 、 退 耕 还 林 ( 还 草 ) . 石 漠 化 治理 等 工程 ,提高 生态 工程 效率 ,实现 生 
态 效 益 最 大 化 。 
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